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ABSTRACT 

We report the results of a globally coordinated photometric campaign to search for transits by the 
P ~ 30 d and P ~ 60 d outer planets of the 3-planet system orbiting the nearby M-dwarf Gl 876. These 
two planets experience strong mutual perturbations, which necessitate use of a dynamical (four-body) 
model to compute transit ephemerides for the system. Our photometric data have been collected from 
published archival sources, as well as from our photometric campaigns that were targeted to specific 
transit predictions. Our analysis indicates that transits by planet "c" (P ~ 30 d) do not currently 
occur, in concordance with the best-fit i = 50° co-planar configuration obtained by dynamical fits 
to the most recent radial velocity data for the system. Transits by planet "b" (P ~ 60 d) are not 
entirely ruled out by our observations, but our data indicate that it is very unlikely that they occur. 
Our experience with the Gl 876 system suggests that a distributed ground-based network of small 
telescopes can be used to search for transits of very low mass M-stars by terrestrial-sized planets. 
Subject headings: stars: planetary systems - stars: individual (Gl 876) - planets and satellites: general 



1. INTRODUCTION 

The Gl 876 planetary system ranks with the most 
remarkable dicoveries that have emerged from the first 
decade of extrasolar planet detections. As described by 
Rivera et al. (2005), Laughlin et al. (2005), Butler et al. 
(2001) Marcy et al. (2001), Delfosse et al. (1998) and 
Marcy et al. (1998), there are a number of reasons why 
the Gl 876 planets are extraordinary. Planet "d" , with its 
1.974 day orbit and Msin(z) = 5.9M ffi is the lowest-mass 
exoplanet known to orbit a nearby main-sequence star 
(Rivera et al. 2005). Planets "c" and "b" are the only 
Jovian-mass companions known to orbit an M-dwarf, and 
their clear participation in 2:1 mean- motion resonance 
gives important dynamical clues to the formation and 
evolution of the system (e.g. Lee and Peale 2002). Fur- 
thermore, the Gl 876 planets are the closest exoplanets 
that have been reliably characterized. Their radial veloc- 
ity amplitudes induce an overall signal-to-noise of nearly 
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100, and they have have been observed for well over a 
decade. 

The detection of Gl 876 d by Rivera et al. (2005) 
underscored the rapid development of the Doppler tech- 
nique. The detection of planets having only a few 
Earth masses is the culmination of advances described 
by Butler et al. (1996), Marcy & Butler (1998), Marcy, 
Cochran, & Mayor (2000), Marcy et al. (2005), and Lo- 
vis et al. (2006). Nevertheless, while radial velocity pre- 
cision has improved to better than 1 m/s, Keplerian or- 
bital fits to Doppler velocities give estimates of M p \ sin(i) 
rather than M p \. In order to ascertain the true mass of 
a planet, one must obtain an independent measure of 
the orbital inclination, i, which is most easily measured 
when the planet is observed to transit across the face of 
the parent star. 

The now-celebrated companion to HD 209458 provided 
the first example of a transiting extrasolar planet, and 
this body has generated an incredible variety of obser- 
vational and theoretical investigations. HD 209458b was 
initially detected with the radial velocity method, and 
was then discovered to transit on the basis of follow-up 
photometric observations at the predicted transit times 
(Henry et al. 1999, 2000; Charbonneau et al. 2000). HD 
209458b's well-sampled light curve, when combined with 
high-precision radial velocity data, has permitted pre- 
cise measurements of the planet's basic attributes (e.g. 
Mazeh et al. 2000, Brown et al. 2001, Laughlin et al. 
2005, Winn & Holman 2005). 

To date, nine additional planets have been found to 
transit their parent stars, and the varied uses to which 
these transits have been put are described in the review 
article by Charbonneau et al. (2006). Our goal in this 
paper is to ascertain whether the two outer planets orbit- 
ing Gl 876 (IL Aqr, GJ 876, HIP 113020) can similarly 
be observed in transit. A positive detection would give a 
further major improvement in our understanding both of 
this specific system and of planetary properties in gen- 
eral. 

Rivera et al. (2005) carried out a photometric search 
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for transits by planet d (P ~ 1.974d), and showed that 
such transits do not occur. To date, photometric searches 
for transits by planet c (P ~ 30d) and b (P ~ 60d) have 
not been reported. Although the a priori geometric tran- 
sit probabilities for Gl 876 c and b are only P tI = 1.4% 
and P t r = 0.9%, respectively, the expected photometric 
transit depths are greater than 10%. Librational motions 
arising from the resonance between the two outer planets 
lead to predicted transit epochs for the outer planets that 
are not spaced evenly in time. (The resonant structure of 
the system is discussed further in Laughlin & Chambers 
2001; Rivera & Lissauer 2001; Kley 2003; and Beauge, 
Michtchenko, & Ferraz-Mello 2005). 

Benedict ct al. (2002) used HST astrometry to derive 
an orbital inclination for Gl 876b of i = 84±6°. Their re- 
sult thus suggests that planet b has a mass equal or close 
to the RV-determined value of Msin(i) = 1.89Mjrjp, 
and that the system is within 1-sigma of being viewed 
edge-on. The implication of the Benedict et al. (2002) 
measurement is that occurence of transits is much more 
likely than the a-priori geometric probability would sug- 
gest. The Benedict et al. (2002) result is in conflict, 
however, with the dynamical models of the system pre- 
sented by Rivera et al. (2005), who find a best fit to 
the radial velocity data when the system (assumed co- 
planar) is inclined by i ~ 50 ± 3° with respect to the 
plane of the sky. Therefore, if transits are detected, then 
the Benedict et al. (2002) results will be confirmed at 
the expense of the conclusions of Rivera et al. (2005). 

Gl 876b and Gl 876c are members of an as-yct un- 
observed class of giant planets with intermediate sur- 
face temperatures. The expected effective temperatures 
of T e ff b = 160k and T c g c = 200K of the planets are 
much lower than the value T e g ~ 1350K measured with 
Spitzer for HD 209458b (Charbonneau et al. 2005), but 
still considerably higher than the T e g ~ 135K measured 
for Jupiter. Determination of their physical properties 
would thus provide a useful link between Jupiter and Sat- 
urn on the one hand, and objects such as HD 209458b 
on the other. Accurate measurement of their radii would 
also give clues to their interior structure (see Charbon- 
neau et al. 2006). 

Finally, the Gl 876 system with its strongly gravita- 
tionally interacting planets and large planet-to-star ra- 
dius ratios, has the potential to display an extraordinar- 
ily informative set of transit light curves. Our model 
light curves herein show what one can expect when in- 
teracting transiting systems are discovered, for example, 
with ongoing radial velocity surveys or with space-based 
missions such as Kepler (Borucki et al. 2003) or COROT 
(Baglin 2003), and allow for a more complete understand- 
ing of the multiplanet dynamics. 

The plan of this paper is as follows. In §2, we de- 
scribe our dynamical model of the system, and illustrate 
some of the photometric properties that the in-transit 
light-curves would be expected to display. In §3, we de- 
scribe an evaluation of the Doppler radial velocity mea- 
surements that have been taken during predicted transit 
intervals. In §4, we evaluate a range of archival and newly 
obtained photometric data sets for evidence of transits. 
In §5, we discuss our results, as well as the current status 
and outlook for the ongoing TransitSearch collaboration, 
and in §6, conclude. 



2. THE DYNAMICAL MODEL 

Gl 876 (M4V) is the fortieth-nearest star to Earth, 
with a Hipparcos-determincd distance of 4.69 pc (Per- 
ryman et al. 1997); Tycho-II and UCAC positions vali- 
date this determination. As discussed in Laughlin et al. 
(2005), current estimates of its mass and radius stand at 
M* = 0.32 M , and R± = 0.3 R Q . We adopt these values 
as fixed in this paper. To date, 155 Doppler velocities 
measured at the Keck telescope and 16 velocities mea- 
sured at the Lick telescope have been published for the 
Gl 876 system (Marcy et al. 2001; Rivera et al. 2005). 
The two planets induce a total velocity half-amplitude 
of the star of nearly 0.25 kms -1 . Furthermore, as dis- 
cussed in Rivera et al. (2005) the system displays a very 
low level of stellar radial velocity "jitter", which further 
aids in obtaining a detailed characterization of the or- 
bits. The system has now been observed for more than 
eighty orbital periods of the middle planet, c. This ex- 
tensive time baseline allows for a much more detailed 
study of planet-planet interactions than can be obtained 
with the other exoplanetary systems known to be in 2:1 
resonance, a list which now includes HD 82943 (Mayor 
et al. 2004, Lee et al. 2006), HD 128311 (Vogt et al. 
2005), and HD 73526 (Tinney et al. 2006). These three 
systems all have inner planet periods of order ten times 
longer than Gl 876 c. Our baseline co-planar i = 90° 
orbital model was obtained by Rivera et al. (2005) from 
a dynamical fit to the Keck radial velocity data. For 
reference, the parameters of this model are reproduced 
in Table l. 13 Additional fits with i < 90° were also ob- 
tained to the Keck data using the algorithm described 
by Rivera et al. (2005). As with the i = 90° baseline 
model, the system is assumed to contain three planets, 
all orbiting the star in the same plane. We varied the 
inclination of the orbital plane to the plane of the sky 
from i = 90° to i = 89.0° in decrements of 0.1°, and fit- 
ted for each of the 13+1 remaining parameters. N-body 
integrations of the fit then generate a sequence of transit 
predictions for the outer two planets. Uncertainties in 
the predicted transit midpoints were generated with the 
bootstrap Monte-Carlo technique described by Rivera et 
al. (2005). For both outer planets, the uncertainty in 
the transit midpoint times was found to vary from ~1 
hr to a few hours. The variation in uncertainty from 
transit to transit is plotted in Figure 1. The predicted 
transit midpoints are best constrained during the epoch 
spanned by the radial velocity observations. Errors in- 
crease as the transit predictions are generated for earlier 
or later times. Light curves were produced with a pho- 
tometric model of the transit which assumes that the 
planets are opaque disks with radii R — 0.93 i?jup for 
c and R = 1.04_Rjup for b. These equatorial radii are 
computed using the models of Bodenheimer, Laughlin, 
& Lin (2003). The specific intensity of the stellar disk is 
modeled with a linear limb darkening law, given as 

I{n)/I{l) = l-v{l-n), (1) 

where I is the specific intensity of the stellar disk, 
v = 0.724 is the corresponding limb darkening coefficient 

13 Note that the elements reported in Table 1 arc osculating 
orbital elements (see, e.g., Murray & Dermott 1999) referenced to 
the epoch JD 2452490. 
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appropriate to V-band observations of an M4V primary 
(per Claret 2000); and fj, = cos(0). 

Figures 2 and 3 show a selection of model light curves 
for planets b and c. For planet b, which has an eccentric 
(e = 0.22), rapidly precessing orbit (dza/dt ~ 44°yr _1 ), 
the shape of the predicted curves depends strongly on 
the osculating periastron angle. For i < 90°, the transits 
are both shorter and deeper when the transits occur near 
periastron. The orbital eccentricity of planet b is more 
nearly circular, with e = 0.025. In this case, there is 
little variation in the predicted light curves as the orbits 
process. The known transiting planets have transit-to- 
transit intervals ATi that are evenly spaced to within the 
current resolution of the observations. HD 209458b, for 
example, has had individual transit midpoints measured 
to an accuracy of several seconds (Brown et al. 2001, 
Wittenmyer et al. 2005), and an average period mea- 
sured to an accuracy of 83 ms (Wittenmyer et al. 2005, 
Winn & Holman 2005). If, however, one is able to mea- 
sure differences in the time intervals between successive 
transits, one can derive additional dynamical informa- 
tion. As discussed by Miralda-Escude (2002), Holman 
& Murray (2005), and Agol et al. (2005), measurements 
of small variations from exact periodicity can be used 
to infer the presence of additional bodies in the system, 
potentially with the masses in the range appropriate to 
terrestrial planets. In the event that either Gl 876b or 
Gl 876c were to transit, one would expect the departures 
from strict periodicity to be large. In fact, successive 
transits of planet c would occur at intervals differing by 
as much as 4 hours! By fitting a dynamical model to such 
transit intervals, one could obtain information regarding 
all of the orbital elements of the planets (including incli- 
nations and nodes). 

3. SEEKING EVIDENCE OF TRANSITS IN THE RADIAL 
VELOCITY DATA SET 

During a planetary transit, the occulting disk of the 
planet passes over portions of the stellar disk contain- 
ing a varying component of line-of-sight velocity arising 
from stellar rotation; this phenomena produces an asym- 
metry in the star's spectral lines known as the Rossiter- 
McLaughlin effect (Rossiter 1924, McLaughlin 1924). 
The effect has been described in detail in the context 
of planetary transits by Ohta, Taruya & Suto (2005), 
and has been analyzed for HD 209458 by Queloz et al. 
(2000), Bundy & Marcy (2000), and Winn et al. (2005), 
and for HD 149026 by Wolf et al. (2006). In fact, the 
transiting planet HD 189733 was initially discovered to 
transit not via photometry, but rather via such in-transit 
radial velocity variations, found by Bouchy et al. (2005). 
For our study of Gl 876, we compared the observation 
epochs of all of the Lick (see Marcy et al. 2001) and 
Keck (see Rivera et al. 2005) radial velocity observa- 
tions to see if any of these observations occurred during 
a predicted transit (assuming i = 90°). For Gl 876, we 
used a rotational period that was determined to be 96.7 
days using rotational modulation (as reported in Rivera 
et al. 2005); we used this value to calculate the am- 
plitude of the Rossitter-McLaughlin effect. Four such 
observations were found, three taken during a predicted 
transit by planet "c" , and one taken during a predicted 
transit by planet "b". Figures 4 and 5 display the rel- 
evant parts of the radial velocity model (both with and 



without the Rossiter-McLaughlin effect), with the obser- 
vations shown. The radial velocity model curves in both 
figures were generated by adding template curves which 
use the stellar and planetary radii and limb darkening 
law given above (and which apply the analytic integrals 
given in Ohta et al. 2005), to our fiducial i = 90° three- 
planet radial velocity model. For the transit epoch in 
the case of planet c (Figure 4), the data is inconsistent 
with a transit. In Figure 5, for the outer planet b, the 
situation is unresolved. 

4. A SEARCH FOR TRANSITS IN PHOTOMETRIC DATA 

SETS 

We have also compared our model light curves to a 
variety of photometric datasets for Gl 876 taken since 
1989. The results of these comparisons are as follows. 

Hipparcos: The Hipparcos Mission (see Perryman et 
al. 1997) produced 67 accepted photometric measure- 
ments for Gl 876 (HIP 113020). These measurements 
have a median magnitude Hp — 10.148, with a = 0.0045. 
There is evidence for a weak periodicity when the data 
is folded at the 96.7 d rotational period of the star (also 
noted in Rivera et al. 2005). In Figure 6, we plot our 
model i = 89.6° light curves for Gl 876 b and c super- 
imposed on the Hipparcos epoch photometry. None of 
the Hipparcos measurements were made near or during 
a predicted transit window. 

Fairborn Observatory: 

A number of long-term photometric monitoring pro- 
grams are being carried out with the automated tele- 
scopes at Fairborn Observatory (e.g., Henry 1999; Eaton, 
Henry, & Fekel 2003) and, as described in detail in Rivera 
et al. (2005), Gl 876 has been on one of the observing 
programs. Figures 7 and 8 show the Fairborn Observa- 
tory data and our corresponding model light curve for 
planets b and c. Data taken during predicted transit 
epochs for planet c show no indication of a transit. The 
Fairborn Observatory data indicate that long-term stel- 
lar variability from Gl 876 is of order 0.05 magnitude, 
which is somewhat smaller than the variation observed 
with Hipparcos. 

Transitsearch and AAVSO: Seagroves et al. (2003) 
describe Transits earch, a cooperative distributed observ- 
ing project involving sub-meter class telescopes world- 
wide. The TransitSearch strategy is to observe known 
planet-bearing stars at the dates and times when tran- 
sits are expected to occur. It is therefore a targeted 
search, which differentiates it from ongoing wide-field 
surveys. We identify observing windows for candidate 
stars through the use of the bootstrap Monte Carlo tech- 
nique described by Laughlin et al. (2005). The observa- 
tional campaigns are prioritized by the a-priori likelihood 
that a particular candidate planet will display transits. 
This likelihood is given by, 



= 0.0045 



1AU \ ( R* + R P i 
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1 + ecos(| — vd) 
l - ^ 

where a is the semi- major axis of the orbit, is the 
radius of the star, R p i is the radius of the planet, w is the 
longitude of periastron, and e is the orbital eccentricity. 

The first TransitSearch dataset for Gl 876 was obtained 
by one of us (Shankland) from 2003 Oct. 27.9583 to 2003 
Oct. 28.0368, with a V-filtered, 16-bit speckle ccd at the 
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primary focus of a 40.4 cm f/5.1 reflector (at air mass 
from sec(z) = 2.7 to 1.7). The photometric time series 
shows a monotonic flux increase of 0.16 +/- 0.04 mag 
from 2003 Oct. 27.9583 to 2003, Oct 27.9972. As shown 
on the left panel of Figure 9, the depth, duration and 
timing of this event were consistent with an egress from 
transit of Gl 876 c, in agreement with our predictions. 
The possibility that a transit had been observed spurred 
us to orchestrate a distributed observational campaign 
during the 2004 season. 

During the 2004 June 24 opportunity, observations 
were attempted by ten observers in Australia. Poor 
weather thwarted all but one observer, and no transit sig- 
nal was seen in the data, which covered about half of the 
1 — a transit window. Eight Australian and Japanese ob- 
servers participated in the following 2004 July 24 planet 
c campaign. They were all clouded out. Weather im- 
proved for the 2004 August 23 planet c opportunity, per- 
mitting observations by five of the nine Australian, South 
African, and German observers. No transit signal was 
seen with photometric depth greater than 1%. 

In 2004 October, a number of AAVSO and Transit- 
Search observers obtained photometry during a closely 
spaced pair of transit windows for planets b and c. Begin- 
ning October 20 11:22 UT and ending October 23 15:27 
UT, a total of 2,795 CCD observations were obtained. In 
addition, 2,981 photometric observations were obtained 
in the days immediately before and after the window to 
set baseline activity and to look for red dwarf flaring. 
Such a large, distributed network of observers allowed 
coverage of the window with very few observing gaps over 
the 3-day period. Gaps in the coverage were smaller than 
the estimated transit period or forecast window, allowing 
us to conclude that transits for planets b and c did not 
occur during the October 2004 campaign. Photometric 
data from the campaign are plotted in Figures 9 and 10. 

Individual observations and uncertainty estimates from 
the October 2004 campaign are available from the 
AAVSO web site, while the other datasets can be found 
at the TransitSearch website 14 . We are continuing to 
obtain photometry of Gl 876, with the goal of gaining 
a better understanding of the photometric variability of 
mature (age exceeding lGyr) red dwarfs. 

5. DISCUSSION 

The TransitSearch photometric network has been in 
operation for three years, and has participating observers 
capable of providing fully global longitude and latitude 
coverage. The network maintains a continuously updated 
catalog 15 of the known census of extrasolar planets. The 
catalog is presently the only available source for (1) a- 
priori transit probabilities and predicted transit depths, 
for the known extrasolar planets, (2) predicted transit 
ephemeris and transit window uncertainties based on or- 
bital fits to published radial velocities, and (3) the results 
of known photometric searches for transits of planet- 
bearing stars. To date, the majority of negative transit 
searches have not been reported in the literature. With 
this paper, we hope to spur a reversal of that trend. 

14 our AAVSO data available at 

http:// www . aavso .org/data/ download / while 
otner collected data can be found at 
http://www.ucolick.org/ %7elaugh/G1876/ c.results.html 

10 available at http://www.ucolick.org/%7elaugh/ 



Due to significant uncertainties in the orbital fits, it is 
generally very difficult to definitively rule out transits by 
planets with P > lOd. The Gl 876 planets form a no- 
table exception to this rule because of their high signal- 
to-noise, r ~ 50, and the extensive radial velocity data 
set allows for a very accurate orbital fit. Campaigns by 
the TransitSearch collaboration have thus far led to tran- 
sits of a number of other shorter-period planet-bearing 
stars being ruled out to a high degree of confidence. The 
planets for which the possibility of transits have been 
significantly discounted by our campaigns include HD 
217107 "b" (P = 7.127d, P tr = 8.0%, see also Vogt et 
al. 2005), HD 168746 "b" (P = 6.403d, P tr = 8.1%, see 
also Pepe et al. 2002), and HD 68988 "b" (P = 6.276d, 
P tr = 9.0%, see also Vogt et al. 2002). 

Furthermore, planets for which usable photometry of 
the parent stars has been obtained by our network during 
the 3 — a transit windows, but for which transits can 
not yet be fully ruled out, include HD 188753Ab (P = 
3.348d, P tr = 11.8%), HD 76700 "b" (P = 3.971d, P tr = 
10.0%), HD 13445 b (P = 15.76d, P tr = 3.4%), HD 74156 
"b" (P = 51.64d, P tr = 3.8%), HD 37605 "b" (P = 
54.23d, P tr = 2.2%), and HD 80606 "b" (P = 111.4d, 
Ptr = 1.7%). Planets for which campaigns are scheduled 
to begin soon include GL 581 "b" (P = 5.366d, P tr = 
3.6%), and HD 99492 "b" (P = 17.04d, P tr = 3.4%). GL 
581 is of particular interest because the primary star is an 
M dwarf, and the planet (if it transits) will be a Neptune- 
mass object M sin(i) = 0.052M Jup (Bonfils et al. 2005). 
The 1.6% transit depth would allow an unprecedented 
physical characterization of a low-mass extrasolar planet. 
To date, our Gl 876 campaign has employed the most 
globally far-reaching distribution, and along with recent 
leaps in fidelity in similar campaigns, clearly underscores 
the utility of the network. 

What is more, Gl 876 b and c constitute a dramatic 
exception to the emerging aphorism that Jupiter-mass 
planets are rarely associated with M dwarfs. Indeed, 
a paucity of giant planets orbiting red dwarfs seems to 
be a natural consequence of the core accretion theory 
of planet formation (Laughlin, Bodcnhcimcr & Adams 
2004). The core-accretion theory predicts, however, that 
Neptune-mass and smaller mass planets will be very com- 
mon around red dwarfs. This stands in contrast to the 
predictions of Boss (2006), who suggests that the gravi- 
tational instability mechanism is the dominant mode of 
giant planet formation. If this is true, then Jupiter-mass 
planets should be just as common around red dwarfs as 
they are around solar-type stars. Efforts to detect tran- 
sits of small planets orbiting small stars are therefore 
likely to gain increasing importance and focus over the 
next several years. 

There are, in fact, several observational indications 
that low-mass planets may be common around red 
dwarfs. The radial velocity surveys have recently re- 
ported the detection of Neptune-mass companions in 
short-period orbits around the red dwarfs GL 581 and 
GL 436. The OGLE team of observers using the mi- 
crolensing method have detected the signature of what 
appears to be a planet with 5.5 Earth Masses orbiting 
a distant red dwarf (Bennet et al. 2006). With these 
points in mind, we can put the results of our Gl 876 
campaign in a broader context. A ground-based (and 
"fiscally viable") photometric network such as Transit- 
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search, especially in collaboration with a dedicated bank 
of telescopes, can capably monitor individual M dwarfs 
to search for transits of terrestrial-sized planets. Owing 
to the intrinsic long-term variability of M dwarfs, candi- 
date planets would have to be identified on the basis of 
a full transit signature in a single time-series. 

A 0.1 M Q M dwarf has R*/Rq ~ 0.1, T cS ~ 2750K, 
and L*/Lq ~5x 10~ 5 ; a habitable planet that receives 
an Earth-equivalent flux from the star thus needs to or- 
bit at a distance of just 0.022 AU, which corresponds to 
an orbital period of 3.85 days. An Earth-like planet with 
such a period would be rotationally synchronized to the 
red dwarf, and in the absence of any significant perturb- 
ing bodies, its orbit would be almost perfectly circular. 
Interestingly however, simulations show that habitability 
is unlikely to be adversely affected by a spin-orbit period 
synchronization. Joshi, Haberlc, & Reynolds (1997) used 
a global climate model to investigate how the Earth's cli- 
mate would respond if the Earth were tidally locked to 
the Sun, and found that Earth remains habitable in this 
configuration - at 1 AU, and perhaps throughout the 
habitable zone. 

The M = 7.5M©, P = 1.9379 d companion to Gl 
876 demonstrates that it is not unreasonable to expect 
terrestrial-mass bodies on short-period orbits of 0.1M Q 
stars. Montgomery & Laughlin (2006) have carried 
out accretion simulations with the Wethcrill-Chambers 
method (Wetherill 1996, Chambers 2001) that model the 
accretion of terrestrial-mass planets in short-period or- 
bits about M dwarfs. These simulations provide further 
evidence that terrestrial-mass planets will commonly 
form in the desired short-period orbits. 

If habitable planets do commonly form in orbit around 
low mass M dwarfs, the chances of detecting and char- 
acterizing them are surprisingly good. The transit of an 
Earth-sized planet will block about 1% of the stellar flux 
of a 0.1M® star. For a planet on a habitable 3.85 day 
orbit, the transit will be relatively brief, ~ 40 minutes. 
The a-priori geometric probability of observing such a 
transit is 2%. 

A 1% photometric dip is readily detectable. Amateur 
astronomers who participate in the Transitscarch collab- 
oration routinely achieve detection thresholds of consid- 
erably better than 1%, as evidenced by confirming de- 
tections of HD 149026b which has a transit depth of just 
0.3%. Indeed, capable amateur observers have demon- 
strated the photometric capability to detect the passage 
of a Mars-sized body in front of an 11th magnitude 0.1 
solar mass red dwarf. 

Inevitably one must ask how many suitable red dwarfs 
are available on the sky as a whole. Even though the 
lowest-mass M stars are the most common type of star, 
they are also exceedingly dim. A dedicated 3m-class tele- 
scope would be required to properly search the magni- 
tude range V=16-18, where observations of the lowest- 
mass red dwarfs begin to be plentiful (and indeed, many 
0.1 M Q stars within 10 parsecs still likely remain to 
be discovered). In the meantime, however, we recom- 
mend that observers obtain high-cadence photometry of 
nearby single stars listed in the RECONS catalog 16 of 
the 100 nearest stellar systems; these include Proxima 

16 maintained by T. Henry and collaborators at Georgia State 
University. 



Centauri, Barnard's Star, Wolf 359, Ross 154, Ross 128, 
DX Cancri, GJ 1061, GJ 54.1, and GJ 83.1. We esti- 
mate that each one has a ~ 1% chance of harboring a 
detectible transiting, potentially habitable planet - read- 
ily detectible by distributed photometric observation. 

6. CONCLUSIONS 

Our study of Gl 876 has produced several interesting 
results. We have demonstrated that transits of planets 
b and c are not currently occurring, and we have thus 
modestly constrained the parameter space available to 
observationally allowed orbital configurations of the Gl 
876 system. The absence of transits is consistent with 
the dynamical study of Rivera et al. (2005) who find 
that a co-planar system inclination i ~ 50° yields the 
best fit to the radial velocity data. Our dynamical anal- 
ysis indicates that when planetary transits in a strongly 
interacting system similar to Gl 876 are eventually ob- 
served, then a tremendous amount of information will be 
obtained. We look forward to the discovery of the first 
such system. 

We have also shown that a globally distributed net- 
work of small-telescope observers can effectively provide 
definitive photometric follow-up of known planet-bearing 
low-mass stars. Our experience with Gl 876 indicates 
that intensive photometric monitoring campaigns of in- 
dividual, nearby M dwarfs constitute a tractable strategy 
for planet detection. The Sun's nearest low-mass stellar 
neighbors should therefore be scrutinized very carefully 
for transits by short-period terrestrial planets. A low- 
budget survey of this type provides a small, but never- 
theless non-negligible chance to discover a truly habitable 
world from the ground, and within a year or two. Such a 
detection would be a remarkable and exciting discovery. 
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Fig. 1. — Uncertainty in the central transit time for planets "b" solid line and "c" dashed line. The vertical dashed lines delineate the 
time span of the Keck RV observations used to generate the co-planar fits that form the basis of our model photometric and RV curves. 
The vertical dotted line indicates the epoch for the fits (JD 2452490.0). 
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Fig. 2. — Model limb-darkened light curves during a predicted transit by planet "b" when it is near periastron (solid curves) and near 
apastron (dashed curves). The curves correspond to different assumed inclinations (i = 90°, 89.9°, 89.8°, 89.5°). The shallowest occur 
for i = 89.5°. 
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Fig. 3. — Model transit light curves for planet "c", when it is near periastron (solid lines), and near apastron (dashed lines). Inclinations, 
i = 90° through i = 89.0° are shown in decrements of 0.1°. Planet "c"'s eccentricity e ~ 0.22 leads to a greater difference in transit depth 
and duration between periastron and apastron. 
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Fig. 4. — Three-planet RV model (dashed curve) with a superimposed model Rossiter curve (solid line) during a predicted transit 
planet "c" near epoch JD 2453301.79. The observed Keck RV data is shown as the black fillcd-in circles with error bars. 
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Fig. 5. — Three-planet RV model (dashed curve) with a superimposed model Rossiter curve (solid line) during a predicted transit by 
planet "b" near epoch JD 2452446.14. Again, the observed Keck RV data is shown as the black filled-in circle with error bars. 
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Fig. 6. — Model i = 89.6° light curves for Gl 876 "b" and "c" superimposed on Hipparcos epoch photometry. None of the Hipparcos 
measurements were made near or during a predicted transit window. 
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Fig. 7. — Model i = 89.3° light curves for Gl 876 "b" and "c" superimposed on Fairborn Observatory photometry. Note that the 
intrinsic long-term variability for Gl 876 over a ~ 1000 day period is of order 0.05 magnitudes. 
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2970.6 2970.8 3271.6 3271.8 3301.6 3301.8 

t (JD - 2450000) 

Fig. 8. — Model light curves for Gl 876 "c" superimposed on Fairborn Observatory photometry taken near three individual predicted 
transit windows. The solid lines correspond to model light curves arising from assumed inclinations running from i = 90.0° to i = 89.0°. 
The light lines correspond to i = 90° model light curves arising from bootstrap fits in which the transit is 3<x early, Icr early, and lcr and 
3<r late. 
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Fig. 9. — Left: Shankland 2003 data for "c", this plate shows dynamical model light curves again for the same decrements of i, taken 
near JD 2452940.4381; relative photometry, which began in astronomical twilight (filled-in circles) Right: Dynamical model light curves 
for Gl 876 due to planet "c" for i = 90.0° through 89.0° in 0.1° decrements, taken near a predicted transit centered at JD 2453301.7879. 
The data were compiled, baselined and shifted to employ a mean magnitude of V=10.19, based on relative photometry obtained by AAVSO 
and TransitSearch (filled-in circles). The light lines correspond to i = 90° model light curves arising from bootstrap fits in which the 
transit is 3<x early, la early, and la late and 3a late. 




Fig. 10. — Dynamical model light curves for Gl 876 due to planet "b" for i = 90.0° through 89.5° in 0.1° decrements, near a predicted 
transit centered at JD 2453300.3806. As with the previous figure, the data were baselined and shifted to employ a mean magnitude of 
V=10.19, based on relative photometry obtained by AAVSO and Transits 'earch(filled-in circles). The light lines correspond to i = 90° 
model light curves arising from bootstrap fits in which the transit is 3a early, la early, and la late and 3a late. 
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TABLE 1 

Co-Planar. Three-Planet i = 90° Jacobi Parameters of the Gl 876 System 



Parameter 


Planet d 


Planet c 


Planet b 




p(d) 

K (ms- 1 ) 
a a (AU) 

« C) 

M (°) 

transit epoch 


5.89 ± 0.54 M 9 
1.93776 ± 0.00007 
6.46 ± 0.59 

0.0208067 ± 0.0000005 
(fixed) 
(fixed) 
309.5 ± 5.1 

JD 2452490.756 ± 0.027 


0.619 ± 0.005 M. 1up 
30.340 ± 0.013 
88.36 ± 0.72 
0.13030 ± 0.00004 
0.2243 ± 0.0013 
198.3 ± 0.9 
308.5 ± 1.4 

JD 2452517.633 ± 0.051 


1.935 ± 0.007 M Jup 
60.940 ± 0.013 
212.60 ± 0.76 
0.20783 ± 0.00003 
0.0249 ± 0.0026 
175.7 ± 6.0 
175.5 ± 6.0 


a Quoted uncertainties 


in planetary masses and semi-major axes do not incorporate the uncertainty in the 


mass of the star 



